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SUGARY 

A  viscous  Interaction  model  has  been  developed  for  the  treat¬ 
ment  of  tangentially  blown  transonic  airfoils.  The  coupling  procedure 
employs  a  small  disturbance  approximation  for  the  nearly  inviscid  outer 
flow  and  interacts  this  region  with  the  boundary  layer,  wall  jets,  and 
propulsive  wake.  Keller's  box  scheme  is  utilized  to  treat  these  viscous 
zones.  On  the  basis  of  parametric  studies  for  laminar  flow,  the  skin 
friction  drag  appears  to  increase  nearly  linearly  with  the  blowing 
coefficient  Cj  for  a  fixed  Mach  number- Incidence  combination,  giving 
a  substantial  scrubbing  drag  penalty  for  "overblowing"  beyond  the  Cj 
necessary  to  achieve  attached  flow.  Increases  in  wave  drag  over  the 
same  range  of  Cj  are  much  smaller.  Associated  supercirculation  gains 
with  overblowing  are  inefficient  from  the  standpoint  of  L/D  due  to  the 
large  values  of  scrubbing  drag.  The  role  of  surface  skin  friction  in 
modifying  the  viscous  cambering  of  the  airfoil  In  the  laminar  case  is 
compared  to  a  turbulent  modification  of  the  coupling  model.  Whereas  the 
laminar  results  show  a  lift  increase  associated  with  the  cambering  action 
of  the  wall  jets,  the  turbulent  solutions  indicate  a  decambering  effect 
resulting  in  decreased  lift.  Also  evident  with  the  turbulent  results  are 
increased  skin  friction  and  wave  drag.  An  explanation  of  these  trends 
in  terms  of  the  Karman  momentum  equation  is  given. 
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NOMENCLATURE 

Slot  blowing  parameters 
Chapman  constant 

Constant  for  pressure  representation  near  nose 

Blowing  coefficient 

Scaled  blowing  coefficient 
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Wave  drag 
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Skin  friction  coefficient 
Slot  height 

Boundary  of  airfoil  surface 

Boundary  of  airfoil  surface,  reduced  stream  function 

Quantity  appearing  in  definitions  of  integral  thicknesses 

Quantity  appearing  in  definitions  of  integral  thicknesses, 
mesh  spacing 

Integral  appearing  in  turbulence  model 
Integral  appearing  in  turbulence  model 

Momentum  flux 

Transonic  similarity  parameter 

Mesh  spacing 

Lift 

Airfoil  chord 
Freestream  Mach  number 

Factor  appearing  in  definition  of c 
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Coefficient  appearing  In  viscous  momentum  equation 
Coefficient  appearing  In  viscous  momentum  equation 

Coordinate  In  direction  of  outward  surface  normal 

Quantity  appearing  In  definition  of  Integral  thicknesses 

Freestream  Reynolds  number 

Coordinate  In  direction  of  parallel  to  surface 

Velocity  component  along  parallels  to  wall,  also  f (n) 

f"(n) 

Cartesian  coordinates 
Angle  of  attack 
Thickness  ratio 
Displacement  thickness 
Wake  displacement 

Eddy  viscosity 
Circulation 

Transformed  vertical  coordinate 
Scaled  slot  height 

Momentum  thickness 

Algebraic  turbulence  model  parameter 
Algebraic  turbulence  model  parameter 
Viscosity 
Density 

Perturbation  potential 
Stream  function 
Vorticity 
Total  vorticity 

Jet  deflection  angle 
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Subscripts 

A  Airfoil 

e  Edge  of  viscous  layer 

u  Upper 

t  Lower 

FF  Far  field 

max  Maximum  value 

SL  Slot 

te  Trailing  edge 

W  Wake 

«»  Freesiream  quantity 

Superscripts 

U  Upper 

I  Lower 

n  Iteration  counter 

Special  Symbols 


<  > 

Average 

[  ] 

Jump 
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1.0  INTRODUCTION 

Achievement  of  separation-free  flow  has  been  an  intrinsic  pacing 
item  in  obtaining  high  lift  characteristics  for  airfoils  and  wings.  At 
low  speeds,  this  capability  is  extremely  important  to  the  Navy  from  the 
viewpoint  of  carrier  take-off  and  landing  performance,  operation  from 
bombed-out  runways  as  well  as  fleet  dispersal  requirements.  From  a 
tactical  viewpoint,  it  is  also  significant  in  high-speed  maneuver  scenarios. 
Associated  with  existing  threats,  combat  in  the  transonic  Mach  number 
regime  will  continue  to  be  important  for  the  indefinite  future.  In  typical 
dogfights,  a  substantial  advantage  is  obtained  with  turn-radius  reduction 
while  maintaining  high  speed  and  excess  thrust.  Accordingly,  aerodynamic 
efficiency  can  play  an  important  role  in  mission  survivability. 

The  attainment  of  high  lift  at  transonic  speeds  is  impacted  by 
shock-induced  separation.  Tangential  blowing  using  engine  air,  ducted 
through  a  chamber, exhausting  gas  along  the  wing  surface  has  been  suggested 
as  a  means  of  coping  with  this  problem,  by  accelerating  the  flow  to  over¬ 
come  the  shock  induced  adverse  pressure  gradients.  In  quantifying  the 
blowing  requirements  to  obtain  separation-free  flow  as  well  as  giving 
information  relevant  to  optimized  aerodynamic  design  of  blown  wings, 
computational  models  can  play  an  important  role.  Although  the  low 
speed  problem  has  received  considerable  attention  since  the  early  work 
of  Spence  in  Ref.  I,  the  viscous-interaction  problem  has  not  been  inten¬ 
sively  analyzed.  However,  recently,  there  have  been  studies  amalgamating 
panel  methods  with  parabolized  Navier  stokes  solvers  to  treat  viscid- 
inviscid  interactions  over  jet-flapped  airfoils  for  the  incompressible 
case  as  exemplified  by  Ref.  2.  To  treat  shock  waves  and  other  nonlinear 
effects  arising  in  the  transonic  regime,  other  procedures  are  required. 

In  this  connection.  Refs.  3-5  give  inviscid  solutions  for  jet-flapped 
airfoils  and  wings  as  well  as  upper  surface  blown  airfoils.  These  solu¬ 
tions  can  be  utilized  as  the  initial  iterate  for  the  nearly  inviscid 
portion  of  the  flow  In  a  Prandtl  weak  interaction  procedure  to  handle 
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attached  cases  appropriate  to  moderate  amounts  of  blowing.  In  Refs.  6 
and  7  the  initial  iterate  of  the  blown  boundary  layer  or  wall  jet  portion 
of  these  flows  is  described.  Therein,  other  theoretical  developments  - 
associated  with  treatment  of  wave  interactions  in  wall  jet  flows  are 
also  provided.  This  work  was  performed  under  the  auspices  of  an  ONR 
Contract  N00014-76-C-0359  during  the  period  January  1976  through  March 
1981  with  the  objective  of  modeling  flows  in  transonic  wall  jets  and 
blown  airfoils. 

This  report  will  describe  the  generalization  of  the  models  given 
in  Refs.  3-7  to  handle  viscous  interactions  over  blown  transonic  airfoils. 
As  such,  it  is  a  summary  of  work  performed  under  ONR  Contract  N00014-81-C- 
0104  during  the  period  16  February  1981  through  15  February  1982.  Con¬ 
tained  in  this  document  is  a  description  of  the  theoretical  formulations 
employed  and  a  discussion  of  results  arising  from  the  computational  flow 
solutions.  The  thrust  will  be  to  treat  the  impact  of  the  coupling  process 
on  the  external  nearly  inviscid  flow  field.  Effects  on  the  airfoil 
characteristics  will  also  be  treated,  including  the  influence  on  the  cir¬ 
culation  of  jumps  in  pressure  and  normal  velocity  in  the  wake  as  well  as 
blown  portions  of  the  airfoil.  Also  included  in  this  analysis  is  the 
attendant  alteration  of  the  far  field  due  to  the  modified  vorticity  in 
the  wake.  With  this  methodology,  the  connection  between  blowing  and 
maintenance  of  separation-free  flow  will  be  discussed  in  the  light  of 
parametric  examples.  Questions  involving  the  impact  of  blowing  on 
achieving  laminar  flow,  minimum  friction  drag,  shock-free  flow  and  opti¬ 
mum  force  characterisitics  will  also  be  addressed. 

Because  of  the  multiple  extrema  inherent  in  the  wall  jet  velocity 
profiles,  difference  rather  than  integral  methods  have  been  utilized 
herein  to  compute  the  viscous  regions  of  the  flow.  Although  integral 
methods  can  represent  substantial  time  savings  in  the  computation,  and 
have  therefore  proven  of  great  value,  their  empiricism,  intrinsic  arbi¬ 
trariness,  and  indifference  to  the  actual  details  of  the  wall  jet  velocity 
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profile  development  can  represent  distinct  disadvantages  In  understand¬ 
ing  blown  flows  with  emphasis  on  tradeoffs  between  separation  control, 
the  different  forms  of  drag, and  the  development  of  high  lift.  In  regard 
to  the  coupling  procedure,  the  method  employed  In  this  analysis  Is  an 
extension  of  those  described  In  Refs.  8-12  to  blowing.  For  solution  of 
the  wall  jet  and  wake,  the  box  scheme  as  described  In  Refs.  13-15  has 
been  adapted  to  solve  the  relevant  Initial  boundary  value  problem  for 
the  boundary  layer  as  It  Is  modified  with  tangential  Injection.  The 
adaptation  to  wall  jets  employs  concepts  that  we  have  developed  In  earlier 
three-dimensional  Incompressible  studies.  An  account  of  these  Is  given 
In  Ref.  16.  In  what  follows,  both  laminar  and  turbulent*  boundary  layers 
are  discussed,  with  emphasis  on  the  former. 

To  treat  the  Inviscid  part  of  the  flow,  the  SLOR  procedures  of 
Ref.  17  which  Is  an  Improvement  over  that  contained  In  Ref.  18  have  been 
employed  In  this  analysis.  The  viscous  interaction  problem  is  solved 
also  by  the  SLOR  procedure  with  modifications  for  the  jump  conditions 
In  the  wake  region.  Details  of  these  are  to  be  found  In  what  follows. 


*The  turbulent  results  were  obtained  under  a  Rockwell  complementary 
IR&D  program  and  are  Included  for  completeness  of  the  discussion. 
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2.0  MODEL  FORMULATION 

A  schematic  of  the  tangentially  blown  transonic  airfoil  to  be^ 
considered  is  shown  in  Figure  la.  An  airfoil  is  shown  with  slots  and 
S_  on  its  upper  and  lower  surfaces  respectively,  which  introduce  air  or 
another  gas  tangential  to  the  airfoil  surface  at  these  locations. 

Regions  @  -(7)  have  special  significance  that  will  be  discussed 
presently.  This  blowing  is  utilized  to  overcome  the  adverse  pressure 
gradient  along  the  surface  introduced  by  shocks  or  the  pressure  recovery 
process  as  the  flow  approaches  the  trailing  edge.  For  laminar  flows, 
and  certain  airfoils,  the  lower  surface  slot  is  usually  required  because 
of  the  enhanced  possibility  of  separation  in  these  cases. 

To  obtain  the  flow  over  this  arrangement,  a  "zonal"  method  is 
utilized  as  contrasted  to  a  full  elliptic  Navier  stokes  procedure.  This 
option  has  been  selected  to  avoid  difficult  grid  generation  issues  and 
lengthy  as  well  as  costly  computational  procedures  to  provide  information 
useful  to  design  and  relevant  to  the  flow  structure. 

For  the  zonal  approach,  a  Prandtl  viscous-inviscid  iteractive 
scheme  is  utilized  to  solve  the  viscous  interaction  singular  perturbation 
problem  for  large  Reynolds  numbers  Re  appropriate  to  practical  applica¬ 
tions.  In  this  limit,  the  viscous  terms  in  the  gasdynamic  equations 
of  motion  are  negligible  except  in  the  boundary  layer.  For  thin  airfoils 
at  small  incidence,  the  Karman-Guderley  (KG)  small  disturbance  equation 
can  be  utilized  to  describe  the  flow  field  in  this  nearly  inviscid  "outer" 
flow  region.  Referring  to  Figure  lb,  the  small  disturbance  assumptions 
can  be  utilized  in  the  usual  way  also  to  transfer  the  flow  tangency 
boundary  conditions  from  the  actual  airfoil  and  wake  surface  to  a  slit 
representing  an  approximate  location  of  the  airfoil's  upper  and  lower 
surfaces  and  its  downstream  projection. 
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Fig.  la  Regions  of  consideration  for  tangentially  blown  airfoil. 


Fig.  lb  Regions  for  interacted  inviscid  boundary  value  problem. 
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2.1  Interacted  Inviscid  Problem 

In  the  viscous-interaction  procedure,  the  inviscid  problem  is 
modified  strictly  through  its  boundary  conditions.  For  the  small  disturb¬ 
ance  framework,  described  here,  the  K6  equation  of  motion  for  the  pertur¬ 
bation  potential  4)  given  by 

[K-(y+1)<|>j^](|)j^x  +  =  0  .  (1) 

valid  inside  the  region  bounded  externally  by  the  contour  ABCD  and 
internally  by  the  slit  POQP  in  Figure  lb  is  solved  subject  to  these 
conditions.  In  (1), 

2/3 

K  =  (1-M„^)/Mj^6  =  small  disturbance  parameter 

6  =  airfoil  thickness  ratio 
=  freestream  Mach  number 

00 

y  =  (M^26)^®y  , 

and  X  and  y  are  Cartesian  coordinates  expressed  in  units  of  the  airfoil 
chord*,  and  subscripts  denote  partial  differentiation. 

With  viscous  interactions,  on  the  basis  of  the  generalization 
of  the  asymptotic  results  in  Ref.  6  to  compressible  flow  or  those 
described  in  Ref.  11,  the  modified  boundary  conditions  can  be  shown 
to  be 


*Slight  corrective  scalings  involving  monomials  in  have  been  introduced 
in  accord  with  empirical  procedures. 
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(fr»(x.0+)  =  F»(x)  +  6  §5^- 


'l>^(x.O-)  =  F' (x)  -  6 


-l^^A 


*{l) 


0<x<l 


dx 


.(2a) 

(2b) 


1/3  3/k 


=  ♦j^(x,0+)  -  4>x(x.O-)  =  -5  (9y  +  6y)e"(x) 


e'(x)  =<(()->=  ^(<fi^(x,0+)  +  (j)~(x,0-)) 


(2c) 

l<x<Xpp  (2d) 

(2e) 


primes  denoting  differentiation  with  respect  to  x  and  where  the  equations 
of  the  airfoil  upper  and  lower  surfaces  respectively  are 


y  =  fifu(x)  -  ax  H  6Fy(x) 


=  6fj^(x)  -  ox  E  6Fj^(x) 


(2f) 

(2g) 


Here.a  denotes  the  angle  of  attack  and  the  wake  deflection  is  given  by 


y  =  6e(x)  .  (2h) 

Also  in  (2),  the  quantities  ^^®  displace¬ 

ment  thickness  on  the  upper,  lower  surfaces  and  wake  respectively,  and 
represents  the  corresponding  momentum  thickness.  These  integral  thicknesses 
are  defined  in  Eqs.  (12)  in  what  follows. 


Additional  boundary  conditions  for  the  problem  (2)  are 
required  in  the  far  field  x*  +  *.  These  are  applied  on  the  finite 

boundary  ABCO  and  are  obtained  under  the  assumption  that  the  viscous 
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far  field  Is  explicitly  of  the  same  form  as  the  inviscid  one,  i.e., 

/’^'(5)itn|z-Cld5  +  f  R(z,5)[d(5)-A]dc 
ttVT  (•'0  *'0 

+  ^R(z.C)  e'(C)dcj 


(3a) 


where 


R(z,C) 


(£+1)*  -  (z+D* 
(C+1)*  +  (z+D* 


(3b) 


A  =  a/5 

and  z  =  x+iv/iTy.  The  quantities  c  and  t  denote  the  camber  and  thick¬ 
ness  functions  of  the  airfoil,  i.e.,  fy(x)  =  c(x)  +  t(x)  and 
and  f^(x)  =  c(x)  -  t(x).  In  addition,  it  is  shown  that 


e'(x)=“VT(  |c^e'(l)  -Ij.g|/2Trx  as  x  -*■  » 


(4a) 


where the  circulation  at  the  trailing  edge  is  determined  numerically 
from 


"te 


[♦] 


x»l 


(4b) 


and  is  a  scaled  blowing  coefficient  such  that  if  Cj  is  the  blowing 
coefficient  defined  as 


Cj  ■  d/qt 
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where  J  is  the  total  momentum  flux  from  the  slots,  q  Is  the  dynamic 
pressure,  l  is  the  chord  then, 

m  113 

S  =  «  <=0  • 

The  expansion  (4a)  Is  used  to  evaluate  the  Infinite  Integral  In  (3a) 
outside  the  computational  domain. 

The  representation  (3a)  Is  used  as  a  Dirichlet  condition 
(p  =  (|)pp  on  ABCO.  Actually,  It  Includes  viscous  Interactions  through 
the  wake  vorticity  which  for  inviscid  flow  according  to  Refs.  1, 
and  5  is  related  to  the  wake  curvature  «e"(x)  by  the  formula 

[\]=fe"(x)  .  (5) 


From  Integration  by  parts,  the  last  Integral  on  the  right  hand  side 
of  (3a)  can  be  related  to  the  wake  vorticity  through  (5).  Actually, 

(3a)  neglects  the  contribution  to  the  vorticity  on  the  airfoil.  This 
has  been  assumed  to  be  a  small  effect  In  our  viscous  Interaction  calcu¬ 
lations.  For  this  reason  and  computational  expediency,  (3a)  was  utilized 
with  the  viscous  effect  occurring  through  updates  of  e'(x)  from  (2e)  in 
the  Prandtl  Iteration  process  Involving  Eqs.  (2). 

A  formula  which  Includes  vorticity  on  the  wing  as  well  as 
In  the  wake  Is  motivated  by  the  picture  that  In  outer  Inviscid  coord¬ 
inates,  the  wake  Is  line  rather  than  a  layer,  and  that  from  Green's 
theorem  as  applied  In  Ref.  3,  the  blown  airfoil-wake  combination  can 
still  be  replaced  by  a  concentrated  potential  vortex  at  the  origin 
x*y*0.  Accordingly,  <j»pp  would  be  given  by 

♦pp  ■  as  X*  +  Ky*  -*•  «> 
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where  " 

00 

-r„  "If  ’ 

•00 

by  Stokes'  theorem,  in  which  r  is  the  circulation  on  a  circuit  at  infinity, 
is  the  vorticity,  is  the  total  vorticity,  and  is  the  jump  in 

potential  across  the  wake  line.  These  ideas  are  consistent  with  the  result 
(for  incompressible  flow)  in  Ref.  20  in  which  the  lift  L  is  given, even 
for  viscous  flow  by  the  Kutta-Joukowski  law  in  which  if  and  U  are 
respectively  the  free  stream  density  and  velocity,  then 

L  =  p„ur„  , 


providing  that  is  the  circulation  on  a  circuit  at  infinity  intersecting 
the  wake  perpendicularly.  Note  that  in  a  viscous  flow,  r  is  circuit 
dependent,  because  of  the  variable  vorticity  along  the  wake.  In  Ref.  21, 
for  a  compressible  flow,  it  is  shown  that 

L  »  -pJXlj  . 


In  the  iterative  cycle  to  be  discussed,  the  first  phase  of  the 
calculation  is  to  solve  the  inviscid  boundary  value  problem.  At  this 
Stage*  6;^  »  and  in  Eqs.  (2a,b,d)  are  set  to  zero  and  the 

coefficient  of  e"(x)  in  (2c)  is  replaced  by  C^/2.  For  this  and  sub¬ 
sequent  Interacted  inviscid  fields,  an  additional  condition  required  to 
make  the  problem  well-posed  is  the  value  of  e'(l).  Also,  (2e)  is  used 
as  indicated  in  Section  3.0  since  the  jet  is  a  free  boundary.  If  t  is 
the  angle  that  the  wake  centerline  or  jet  leaves  the  trailing  edge,  then 
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€'(1)  *  -(a+T)/6  .  (6) 

For  a  jet  flap,  i.e.,  the  jet  Is  at  the  trailing  edge,  and  t  is  specified. 
For  tangential  blowing,  the  inviscid  condition  is  applied  that  the  wakS* 
centerline  for  a  wedge  shaped  trailing  edge  leaves  tangent  to  the  higher 
stagnation  pressure  surface, which  for  the  cases  considered, is  on  the 
upper  side.  Accordingly, 

t  =  -1806  f*  (l)/iT  (7) 

is  assumed.  In  actuality,  the  lift  and  some  flow  details  could  be 
affected  by  this  choice  and  more  accurate  "multi-deck"  viscous  local 
solutions  generalizing  those  of  Stewartson  and  other  workers  (see 
Ref.  22}  would  appear  warranted  to  refine  (7). 

The  foregoing  describes  essential  features  of  the  calculation 
of  the  viscous  corrections  to  the  primarily  inviscid  portion  of  the 
blown  flow.  To  complete  this  part  of  the  formulation,  viscous  correc¬ 
tions  to  the  surface  pressures  are  indicated.  These  are  implemented 
at  the  conclusion  of  the  viscid-inviscid  iteration  sequence  and  have 
the  forms 

[<|)jj(x,0+)]  =  (i)j^(x,0++)  -  (>j^(x,0+)  =  (8a) 

[(D'^(x,0-)]  -  ♦^(x,0-)  -  <Dj^(x,0-)  =6^'*  F;^(x)  (8b) 

for  0  <  X  <  1 

where  -h-  and  —  denote  the  edges  of  the  boundary  layer  or  wall  jet  over 
the  upper  and  lower  surfaces  respectively,  and  +  and  -  indicate  the 
upper  and  lower  surfaces  respectively.  At  x  -  0  in  the  computations, 
the  quantities  zero. 
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2.2  Interacted  Viscous  Problem 

To  solve  the  Inviscld  problem,  the  Integral  thicknesses, 

I.e.,  6*'s  and  6's  appearing  In  Eqs.  (2)  must  be  determined  from  the  - 
viscous  region  of  the  flow.  For  this  purpose,  the  adaptation  of  the 
box  scheme  to  tangential  blowing  as  described  In  Ref.  6  has  been  used. 

To  achieve  the  necessary  clustering,  a  stretched  vertical  coordinate  n 
Is  Introduced.  If  (s,n)  represents  a  coordinate  system  erected  respect¬ 
ively  along  parallels  and  normals  to  the  airfoil,  then, 

n  =  p(Ug/pgUgS^)*n  (9) 


(n  Is  the  direction  of  the  outward  normal)  where  p  Is  the  density,  u  Is 
the  velocity  component  parallel  to  the  surface,  subscript  e  denotes  the 
edge  of  the  boundary  layer,  wall  jet  or  wake,  y  Is  the  viscosity,  and 
refers  to  the  distance  along  the  surface  from  the  stagnation  point.  For 
thin  airfoils,  n  «  |y|,  and  x  away  from  the  nose  region  of  a  blunt 
nosed  airfoil.  These  approximations  will  also  be  assumed  valid  on  the 
wake.  If  the  stream  function  ip  Is  represented  as 

V  =  (PeVe^^*  f(n.x)  (10) 


then  the  velocity  ratio  u/Ug  *  =  f'(Ti,x).  This  leads  to  the  transformed 

momentum  equation 


(bf-y 

where 

=  Pe/P  “  ’  "•  =  (t)  /Ci+2/(Y-l)M*] 


(11a) 


(11b) 
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nij  =  (1/2)  (1+1112)  (lid) 


b  E  C(l+et) 

(lie) 

C  =  (pu)/pgUg 

(Ilf) 

Here,  is  an  edc(y  viscosity  which  will  be  discussed  later.  For  laminar 

flow,  b  -  1.  In  the  solution  procedure  and  in  (11b),  an  adiabatic  wall 
and  constant  stagnation  enthalpy  in  the  blown  boundary  layer  and  wake 
has  been  assumed.  Furthermore,  a  Chapman  gas  approximation  in  which 

C  -  1  has  also  been  made. 

From  this  formulation,  the  integral  thickness  appearing  in 

Eqs.  (2)  are  given  by 

*  'jJ 

[“sdn 

0 

(12a) 

6*(^)  =  qj 

*  i 

gdn 

Loo 

r“ 

(12b) 

-Oj 

gdn 

Loo 

(12c) 

e/“)  .  qJ 

[“hdn 

(12d) 

.  QJ 

fO 

hdn 

Loo 

(12e) 

Ou  *  Q 

J 

r hdn 

Loo 

(12f) 
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where 


g  H  c-f  (12h) 

h  =  f(l-f')  (12i) 

Re  =  freestream  Reynolds  number 

For  numerical  purposes,  the  infinite  upper  limts  appearing 
in  Eqs.  (11)  and  what  follows  are  replaced  by  ±n„,  where  n^»l  represents 
edge  of  the  boundary  layer,  wall  jet,  or  wake. 

In  accord  with  the  discussion  in  Ref.  6,  the  parabolic  initial 
boundary  value  problem  for  Eq.  (11a)  involves  the  Following  conditions: 


f(0,x)  =  0 

(13a) 

f'(0,x)  =  0 

(13b) 

f(<»,0)  =  1 

(13c) 

f'(-«,0)  =  1 

(13d) 

Over  the  upper  surface  and  lower  surface  of  the  airfoil,  the  boundary 
layers  in  regions  (D  +  (2)  in  Figure  la  are  solved  subject  to  (13a, b) 
as  well  as  (13c)  and  (13d)  respectively.  Initial  conditions  for  regions 
(!)-<■  (2)  are  obtained  by  specializing  Eqs.  (11)  to  x  =  0  and  solving 
the  resulting  two  point  boundary  value  problem.  In  this  process,  the 
asymptotic  solutions  in  the  neighborhood  of  a  finite  curvature  airfoil 
nose  in  the  K-G  small  disturbance  theory  as  well  as  those  for  viscous 
stagnation  flow  are  used  to  obtain  the  various  coefficients  appearing 
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in  (11a)  given  by  Eqs.  (11b)  -  (lid).  Further  details  of  this  procedure 
are  given  in  Section  3.0. 

The  box  scheme  integration  as  described  in  Section  3.0  is  con¬ 
tinued  until  the  slot  locations  x  =  reached.  There¬ 

upon,  new  initial  conditions  are  applied.  For  the  NACA  0012  laminar 
cases  considered,  a  slot  and  tangential  injection  was  required  on  the 
lower  surface  to  avoid  separation  there  and  thereby  permit  a  fully  con¬ 
verged  Prandtl  iteration  sequence.  As  outlined  in  Ref.  6,  the  initial 
condition  of  the  slots  for  which  "  ^SL  assumed  for 

the  parametric  studies  employed  the  marched  boundary  layer  from  upstream 
(x  <  Xj^)  for  n  >  Tij|»  and  an  injection  profile  which  was  assumed  parabolic*, 
i.e. 


.  0<lnl<lndl  ,  (I4a) 

SL 

where  a  translation  of  axes  is  made  such  that  n  =  0  is  the  position  of  the 
airfoil  surface  downstream  of  the  slot,  u  and  SL  refer  to  upper  and  lower 
surfaces  respectively,  and  the  slot  location  is  given  by  n  =  d.  In  (14), 
the  constants  A^  andn^j  are  related  to  the  blowing  coefficient  as  follows: 

A 

Cj-  (n5)(A|.Aj)(xs|_/Re)*/15  (14b) 

.  (16/15)(<i/4)(p/p,)(f;<;>%f;jJ>') 


where 

(14c) 
(14d) 

*More  general  distributions  could  be  investigated  with  our  computational 
model. 


•^d  *  c  “  transformed  slot  height 


■‘f 


p  *  d  I  pdy  >  mean  density  across  slot 
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and  the  subscript  max  refers  to  the  maximum  value.  The  quantity  p/p^ 
above  can  be  used  to  obtain  a  temperature  effect  in  modifying  Cj.  In 
the  calculations  reported  in  Section  4.0,  the  blowing  parameters  of  the” 
upper  and  lower  slots  were  made  to  be  the  same.  Furthermore,  A^, 
values  were  sought  giving  attached  flow  over  the  entire  airfoil  surface. 

However,  when  separation  did  occur  for  an  assumed  value  of  the  they 

£ 

were  increased  until  fully  attached  flow  conditions  were  achieved. 

Results  for  blowing  inadequate  to  maintain  completely  attached  flow 
have  been  used  to  correlate  the  separation  point  x^^^  against  Cj  in 
Refs.  6  and  7  as  well  as  in  Section  4.0. 

To  treat  the  flow  between  the  slot  locations  x  =  x^^^  and  the 
trailing  edge,  i.e..  Regions  Q)  in  Figure  la,  the  initial  distributions 
f"  and  f  in  our  scheme  are  obtained  by  analytic  differentiation  and  inte¬ 
gration  to  supply  the  portion  of  the  initial  data  on  the  line  x  =  x^j^ 
lying  in  the  range  0< |ri|<|n^^l.  As  indicated  previously,  the  part  of  the 
initial  data  on  x  =  x^j^  lying  in  the  range  |njj|<lnl < ln„|  is  obtained 
from  box  scheme  values  at  x  =  x^^^  associated  with  the  integration  of 
(11a)  from  x  =  0  to  x^^^.  For,  Regions  Q)  ,  the  integration  process 
continues  until  the  trailing  edge  x  =  1  is  reached.  There,  the  integ¬ 
rated  solution  of  the  upper  surface  is  combined  with  that  coming  from  the 
lower  side  to  provide  initial  data  on  the  interval  for  downstream 

integration  of  the  wake-jet  (Region  ®  in  Figure  la).  The  integration 
continues  until  the  far  field  boundary  BC  in  Figure  lb  is  reached,  i.e., 

X  =  Xpp.  The  boundary  conditions  used  in  Region  ®  are  (13a, c,d).  Our  numerical 
approach  to  satisfy  them  is  described  in  Section  3.0. 

Algebraic  Eddy  Viscosity  Model 

For  the  turbulent  results  described  in  Section  4.0,  the  following 

23 

two  layer  edc(y  viscosity  model  of  Patankar  and  Spalding  as  well  as 
24 

Ramaprian  was  us«;d. 
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First  Layer 

cj,  =  (K^V^lf"(Ti)l)I/c2  0<y<y* 

Second  Laver 

ej,  -  a*VxR61f''(n)l)I|/c^ 

where 


Ij  =  Kni) 

y*  = 

and  is  the  first  n  (beginning  at  ri=  and  proceeding  to  n  =  0)  where 


ll-f'(ni)l  <0.01 


The  quantities  k  and  X  in  (15)  are  empirical  constants. 

In  Regions  O,  we  used  typical  wall  jet  values  of  X  =  0.125 
and  <  =  0.435.  For  0<x  <  0.005,  a  laminar  model  was  used.  For  other 
regions,  we  used  X  -  0.1  and  k  -  0.4.  More  sophisticated  models  have 
been  proposed  involving  intermittency  and  damping  corrections  as  well 
as  a  special  procedure  for  wakes.  The  results  in  Section  4.0  do  not 
contain  these  refinements.  Instead,  Eqs.  (15)  are  used  to  address 
some  of  the  more  global  trends  and  issues  that  arise  in  comparing  the 
turbulent  case  against  the  laminar  one,  particularly,  regarding  the 
effect  of  blowing  on  aerodynamic  characteristics  such  as  lift  and  drag. 


(15a) 


(15b) 
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Our  computational  model  however  Is  general  enough  to  assimilate  other 
turbulence  simulations. 

We  note  that  In  the  Implementation  of  (15)  in  the  wake  regionr 
two  Intervals  of  validity  of  Eq.  (15b)  are  used,  I.e.,  y..<y<®  and 
-'»<y<yjj^.  In  a  central  region,  (15a)  Is  used  for  which 
Other  aspects  of  the  application  of  (15)  Involve  assuming  that 
c'(tj)  =  0  for  X  stations  close  to  x  =  Xj|^  and  x  =  1.  This  measure  appears 
necessary  to  stabilize  the  box  scheme.  Additionally,  n„  must  be  made 
larger  In  the  turbulent  case  to  handle  the  more  rapid  growth  of  the 
turbulent  layer  than  the  laminar  one. 
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3.0  NUMERICAL  PROCEDURES 


Inviscid  (0^*^  Iterate)  _ 

The  numerical  method  for  solving  the  interacted  inviscid  problem 

comprised  by  Eqs.  (1)  -  (3)  is  the  standard  SLOR  of  Jameson^^  which 
18 

includes  Murman's  fully  conservative  relaxation  (FCR)  technique.  In 
the  subsonic  region,  centered  differences  are  employed,  while  in  the 
supersonic  zone, upwind  differencing  is  necessary  for  stability.  The 
boundary  conditions  on  the  airfoil  and  the  far  field  are  incorporated 
in  the  usual  way.  For  the  wake  region,  Eq.  (5)  is  integrated  to  yield 

[<l>]  =  r^g  +  C^l  e'(x)  +  (a+T)/6|  /2  (5") 

where  (4b)  and  (6)  were  used. 

In  our  procedure,  Eq.  (S')  is  inserted  into  the  discretization 
0^  3t  the  line  above  y  =  0  by  replacing  (|».  .  .with 

jrjr  I 

‘*’i,j-l  ^  ^e  ^  (a+T)/6}  /2  ,  (16a) 

and  at  the  line  below  y  =  0  by  replacing  41^  with 

**'i,j+l  ■  ^e  ■  S  {  ^  /2  .  (16b) 

for  l<x^<Xpp. 

Further  details  are  given  in  Ref.  4.  After  each  sweep  of  the  grid  network, 
e'  is  updated  using  Equation  (2e),  i.e., 

(nl  .  /  \ 

e'  (x^)  =  j  '♦i,j-l)) 

where  y^  is  the  first  mesh  point  below  the  line  y  «  0  and  the  superscript 
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n  denotes  the  relaxation  sweep  counter.  Unfortunately,  because  of  the 
singularity  at  the  trailing  edge,  Eq.  (17)  cannot  be  applied  for  the 
first  two  mesh  points  to  the  right  of  x  ®  1.  Here,  e'{x^)  and  e'Cxg)  are 
determined  using  an  Interpolation  procedure  and  Eq.  (6),  I.e.,  In  which 
e'j,  and  are  Interpolated  using  eg=  e'  (1)  =  -(a+T)/6  as  well  as  e^, 
e5i,  ec.  and  eU,  which  are  determined  from  (17).  Convergence  occurs  when 

4  b  0  _c 

the  maximum  residuals  for  (p,  e'  and  I^g  are  all  less  than  10  . 

Viscous  Zones 

The  boundary  layer  equations  and  boundary  conditions,  (11)  and 
(13)  are  discretized  using  Keller's  box  scheme  (13,  14).  First,  these 
equations  are  re-written  as  a  first  order  system  by  introducing  the  new 
Independent  variables  u(ri,x)  and  v(ri,x),  giving 

f  =  u  (18a) 

u'  =  V  (18b) 

(bv)'  +  m^fv  4ni2(c-u^)  =  x  (u  |^  -  v  (18c) 

Using  Keller's  box  scheme  on  a  nonuniform  grid,  (18a,  18b)  are 
discretized  as 


(19a) 


(19b) 
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where  we  have  averaged  about  the  point  Using  the  box  centered 

at  (18c)  Is  discretized  as: 


(bv)'l  -  (bv)"  , 


/f"  - 

rj-i 

j-4  /  + 

n-l^n 

h-i 

kn 

•^n  J  "-1 

(19c) 


where 


Vl 


(bv)r[^_(^ 


n-1 


"  Vi  (“j-K-i  ''j-i^j-i)/'^n  • 

where  h.  =  i*  and 

w  J  w  * 

"r.  *  ’‘n  -  “n-r  “">*  '"-»  *  *<^3  *  'j-l’  ■  '“=■ 

The  boundary  conditions  become 


(19d) 


fJJ  =  u"  =  0 

0  0 

(19e) 

“J  “  1 

(19f) 

since  the  parameters  m^  and  m2 
vicinity  of  the  singularity  at  the  nose 

are  difficult  to  compute  In  the 
of  the  airfoil  (x“0).  It  Is 
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necessary  to  develop  the  following  local  solution*  applicable  to  finite 
leading  edge  curvature  noses. 


2/3  _  2/3  2/3 

u  1-6  Cf»/(V +  6C«  )  ,  0<x*<Xr 


assuming  U  =  1,  A  =  0(1),  and  where 


(20a) 


X*  =  (x/ 1.48456) 2 


(20b) 


for  the  NACA  OOXX  series  airfoils.  The  quantity  Cq  is  obtained  by  fitting 
with  the  numerical  solution,  i.e.. 


Mco"*Cp(xc_  )  X(.; 


_£L 

2 


(20c) 


where  Xq  is  the  first  mesh  point  to  the  right  of  x=0  and  Cp(Xpj)  is 

the  corresponding  pressure  coefficient  (computed  numerically)  at  this 
mesh  point.  For  x  >  x^.  a  cubic  spline  fit  is  sufficient  for  interpola¬ 
tion  of  Cp  data  where 

3/h 

for  X*  >  x-  .  Using  (20)  or  (21),  m,  and  m,  may  be  easily  computed  from 

Si  ^ 

(11c)  and  (lid). 


♦Note  that  for  A"^  =  o(l),  according  to  Ref.  25,  (fr,  (sgn  y)x*"^^^  and 


Cq  must  change  sign  accordingly. 


22 


I 


Rockwell  International 

Sci«nc«  C«ntMr 


SC5296.1FR 

The  discretized  system  of  equations  (19)  may  be  solved  very  efficient¬ 
ly  using  Newton's  method  because  the  Jacobian  matrix  has  a  block  tri-diagonal 
structure.  Furtermore,  the  box  scheme  is  second  order  in  a  variable  gcid 
and  can  handle  discontinuities  in  the  coefficients  m^  and  m2>  which  are 
very  important  at  x^j^,  the  shock  location,  and  the  trailing  edge.  See 
Refs.  13  and  14  for  the  details. 

As  indicated  in  Ref.  26,  if  the  initial  condition  is  not  the 
natural  one,  then  physical  disturbances  and/or  spatial  oscillations  may 
occur  in  the  computed  solution.  At  the  slot  location,  it  is  necessary  to 
damp  out  these  oscillations  by  using  a  fine  grid  (Ax  =  0.001)  and  averaging 
over  the  previous  solution  (Xj^_j)  for  the  first  five  mesh  points  to  the 
right  of  x^^.  This  procedure  works  well,  but  has  the  effect  of  decreasing 
the  order  of  the  local  truncation  error  in  the  discretized  system  from 

second  to  first  order.  Since,  however,  the  mesh  has,  also  been  substan¬ 
tially  refined  no  loss  of  accuracy  should  result. 

In  the  wake,  a  three  point  boundary  value  problem  can  be  solved 

at  each  x  station  to  satisfy  (13a, c,d).  However,  for  coding  expediency, 

as  outlined  in  Ref.  6,  we  used  an  iterative  method  related  to  a  suggestion 

27 

of  Cebeci  and  other  iterative  schemes  we  considered  in  which 
f (-njj* x)  =  y(x)  is  iterated  on  y  until  (13a)  is  satisfied.  Here,  y(x)  is 
determined  by  Newton's  method  to  satisfy  (13a).  The  process  usually 
converges  in  one  to  three  iterations. 

More  sophisticated  treatments  are  possible  in  which  Goldstein- 
like  layers  are  interposed  near  the  wake  centerline  to  handle  the  discontin¬ 
uity  in  boundary  conditions  at  y-0  in  proceeding  from  the  left  to  right 
across  x=  1, and  the  associated  secondary  layer  near  the  wake  centerline. 

An  account  of  one  of  them  for  synnetric  wakes  is  given  in  Ref.  28. 

Interacted  Inviscid  Region 

The  technique  for  solving  the  interacted  inviscid  region  is 
similar  to  the  SLOR  method  previously  described.  The  airfoil  conditions 
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(2)  are  handled  similarly.  However,  on  the  wake,  the  jump  conditions  (2c, d) 
are  Incorporated  and  require  special  attention.  No  mesh  points  appear  on  the 
line  y  =  0.  The  first  points  above,  and  below  the  axis  are  at  ±h/2  and- 
the  second  points  are  at  ±3h/2.  By  Taylor's  series  and  (2c)  and  (2d) 


<l>(x,|)  -  (J>(x,^)  =  [4»(x,0)]  +  I  U~(x,0)]  + 


1/3  3/4  rX  * 

I^e  -  ^  Moo  (Qw  "  ^w) 


y' 

dx 


(22) 


J] _ w 

26  dx 


where  the  Integral  may  be  evaluated  numerically.  Equation  (22)  may  now 
be  inserted  Into  the  discretized  form  of  (|)~~  using  the  same  "trick"  given 
by  (16a)  and  (16b). 

Equation  (2e)  is  discretized  above  and  below  the  y-ax1s,  and  the 
resulting  expression  replaces  Eq.  (6)  for  x>  x^.  For  Xj  and  X2,  the 
previous  interpolation  technique  must  be  used  because  of  Instabilities 
created  by  the  singularity  at  the  trailing  edge. 

The  algorithm  described  above  works  quite  well  for  laminar 
problems  and  no  underrelaxation  Is  ever  necessary,  in  contrast  to  the 
Integral  methods  described  by  Lock  In  Ref.  11.  As  would  be  expected,  the 
Integral  methods  for  solving  the  boundary  layer  problem  are  faster;  however, 
the  boundary  layer  profiles  are  much  more  complicated  (multiple  ext>^ema) 
when  blowing  Is  Introduced,  and,  consequently,  we  believe  that  Integrating 
Eq.  (14)  using  the  box  scheme  Is  a  more  accurate  approach.  The  only  dif¬ 
ficulties  In  our  method  sometimes  occur  at  the  slot  location,  trailing 
edge  or  at  the  shock  locations  at  high  Mach  numbers.  The  remedy  Is  to  use 
a  finer  grid  In  these  areas  and  to  smooth  the  cusped  velocity  profiles 
(f)  at  the  slots  and  trailing  edge. 
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Convergence  of  the  overall  problem  Is  defined  to  occur  when  the 
C|_'s  from  two  complete  viscous-inviscid  cycles  converge  to  0.001.  A  flow 
diagram  of  our  algorithm  appears  In  Figure  2.  _ 


t 
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F1g.  2  Flow  diagram  of  algorithm. 
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4.0  RESULTS 

To  illustrate  the  effect  of  viscous  interactions  on  the  perfocm- 
ance  of  tangentially  blown  airfoils,  the  coupling  procedure  described  in 
previous  sections  was  applied  to  a  NACA  0012  airfoil  primarily  within 
laminar  flow  environments.  In  this  stu^y,  the  influence  of  such  param¬ 
eters  as  Mach  number,  incidence,  blowing  coefficient,  Reynolds  number 
and  slot  position  on  the  aerodynamic  characteristics  was  investigated. 

A  solution  library  including  approximately  25  cases  has  been 
developed  which  has  rather  extensive  information  regarding  the  flow 
fields  for  these  varying  conditons.  From  it,  typical  results  have  been 
selected  and  discussed  herein  to  highlight  trends  and  issues  of  engineer¬ 
ing  significance.  It  will  be  clear  that  they  are  also  relevant  to  more 
fundamental  concerns  regarding  flows  over  blown  wings.  Some  of  these  are: 

•  The  structure  of  jet-wakes,  particularly  their  effective 
velocity  and  impact  on  the  global  supercirculation  and  lift 

•  The  genesis  and  control  of  wave  and  friction  drag  in  blown 
transonic  flows 

•  Growth  properties  and  momentum  budgets  of  wakes  from  pro¬ 
pulsive  wings,  invariances  associated  with  these  budgets 
and  deviations  from  self  similitude  for  coflowing,  jets  as 
compared  to  free  and  self-propelled  jets 

•  Blown  Kutta  condition 

•  Stability  and  transition  of  wall  jets  and  tangentially 
blown  flows  over  airfoils  and  wings 
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a  The  structure  of  the  nearly  inviscid  far  field  and  its 
asymptotic  matching  to  the  near  field 

a  Multiple  decks  near  the  jet  exit  and  trailing  edge,  break¬ 
down  of  boundary  layer  approximations  and  wave  interactions 
inside  wall  jets  and  mixing  layers. 

Because  separation  phenomena  were  not  considered  in  the  analysis, 
the  results  for  the  unblown  case  will  not  be  included  in  this  discussion. 
For  this  reason,  and  as  suggested  previously,  it  was  also  necessary  to 
introduce  blowing  on  the  lower  surface  to  overcome  the  adverse  pressure 
gradients  typically  present  for  this  particular  airfoil,  even  at  zero  inci¬ 
dence.  Configurations  with  blowing  on  the  lower  surface  as  well  as  the 
upper  one  relate  to  transonic  applications  of  augmenter  wings  and  other 
schemes  for  high  lift  development  at  these  Mach  numbers. 

A  question  of  some  importance  in  the  utilization  of  the  coupling 
scheme  indicated  in  Figure  2  is  the  number  of  viscous-inviscid  interations 
required  for  convergence  of  the  coupling.  In  Figure  3,  chordwise  pressure 
distributions  for  a  typical  laminar  case  are  shown  for  the  NACA  0012  air¬ 
foil  at  M^  =  0.75,  a  =  2“,  and  a  blowing  coefficient  C.  =  0.04  with  the 

00  J 

upper  and  lower  slots  at  19%  chord.  As  indicated  in  Figure  2,  the  initial 
iterate  is  the  inviscid  blown  solution.  This  gives  rise  to  the  pressures 
shown  in  Figure  3.  Subsequent  couplings  move  the  terminating  shock  for 
the  supersonic  region  rearward.  The  results  show  that  the  convergence  of 
the  viscous-inviscid  cycle  is  extremely  rapid,  requiring  only  two  itera¬ 
tions  in  this  case.  Two  to  three  iterations  was  typical  for  all  of  the 
approximately  25  cases  studied.  In  this  connection,  it  is  interesting 
and  important  to  note  that  the  introduction  of  blowing  provides  a  new 
opportunity  to  stucty  weak  interactions  not  generally  accessible  to  dif¬ 
ference  methods  for  the  transonic  unblown  case,  unless  inverse  procedures 
are  used  and  other  special  measures  such  as  backward  marching  and  degen¬ 
erate  forms  of  the  boundary  layer  equations  are  employed  to  treat  the 
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attendant  separated  regions.  The  significance  of  this  aspect  is  that  it 
may  allow  investigation  of  the  treatment  of  the  trailing  edge-near  wake 
region  in  an  unseparated  context  and  therefore  enable  us  to  understands^ 
its  matching  to  the  external  flow  more  easily. 

Another  feature  of  the  laminar  distributions  shown  in  Figure  3 

is  the  fact  that  they  exemplify  a  case  in  which  the  viscous  lifting 

pressures  are  actually  higher  than  the  "inviscid"  values.  This  somewhat 

surprising  result  which  differs  as  will  be  seen  from  the  corresponding 

turbulent  flow  contrasts  with  the  unblown  case  in  which  viscosity  reduces 

the  lift.  Interpretation  of  this  finding  relates  partially  to  the  viscous 

character  of  the  inviscid  blowing  model.  This  can  be  appreciated  from 

the  asymptotic  thin  jet  analysis  given  in  Ref.  5  which  is  a  generalization 

and  asymptotic  formalization  of  that  given  by  Spence  in  Ref.  1.  In  this 

connection,  it  appears  questionable  whether  the  latter  is  the  infinite 

Reynolds  number  limit  of  the  weak  interaction  theory  formulated  in  this 

report.  Another  more  important  factor  relates  to  the  change  in  sign  of 

the  wall  jet  displacement  and  momentum  thicknesses  and  their  streamwise 

slopes  associated  with  blowing  as  compared  to  a  conventional  boundary 

layer.  According  to  the  transpiration  boundary  conditions  (2a, b),  this  can 

result  in  an  effective  cambering  and  incidence  change  induced  by  the  wall 

jets  and  jet  wakes  as  contrasted  to  the  decambering  action  of  conventional 

boundary  layers  for  the  unblown  case.  A  qualitative  discussion  of  this 

mechanism  is  given  in  Ref.  20.  Along  with  the  cambering  tendency,  the 

negative  displacement  thicknesses  result  in  thinning  of  the  profile  as 

contrasted  to  thickening  occurring  with  ordinary  boundary  layers.  The 

former  is  illustrated  in  Figure  4,  where  displacement  and  momentum 
★ 

thicknesses  6^  and  0^  over  the  upper  su.'face  of  the  NACA  0012  airfoil 

are  shown  along  with  the  local  skin  friction  coefficient  and  peak  velocity 

variations  for  the  M  =  0.75,  a  =  1°,  C,.  =  0.055  case.  At  the  slot 

location,  x  »  0.19,  the  quantities  jump,  related  to  the  mass  flow  from 

★ 

the  slot.  Downstream  of  this  location,  the  displacement  thickness  6^ 


30 


SC82-17414 


Rockwell  International 

Sctonc*  C«ntor 


SC5296.1FR 


I  I 


AlllNVnO 


31 


Variation  of  integral  thicknesses,  peak  velocity  (f'max) 
skin  friction  coefficient  Cj  for  Moo=  0.75,  a=  1°,  C,  »  0.055 
on  upper  surface  of  blown  NACA  0012  airfoil.  “ 


obtained  from  Eqs.  (12)  is  negative,  related  to  the  excess  velocity  wall 

jet  which  contrasts  to  the  defect  boundary  layer  encountered  in  the  unblown 

★ 

case.  In  addition,  negative  <S^'(x)  variations  are  apparent  between  the 
slot  and  trailing  edge  relevant  to  the  aforementioned  sink-like  action. 
Associated  with  this  behavior  of  the  integral  thicknesses  is  the  stream- 
wise  evolution  of  the  velocity  profiles  f'(n)  depicted  in  Figure  5.  Here, 
the  plus  and  minus  signs  refer  to  conditions  immediately  downstream  and 
upstream  of  the  slot,  respectively.  Analogous  development  is  shown  in 
Figure  6  for  the  trailing  edge  and  wake  (x  >  1)  profiles.  Although  not 
attempted  in  this  investigation,  this  information  can  be  used  to  perform 
a  study  of  the  effect  of  the  viscous  interaction  on  transonic  thrust 
recovery,  and  has  a  bearing  also  on  the  development  of  the  wake  vorticity, 
curvature,  and  circulation  as  they  affect  lift.  In  Ref.  6,  the  asymptotic 
development  of  the  wake  profile  to  a  nearly  universal  bell  shape  profile 
is  discussed.  Similarities  and  differences  of  the  present  application 
to  flat  plate  wakes  as  well  as  free  and  self-propelled  jets  are  described 
therein. 

In  Figure  7,  the  effect  of  blowing  coefficient  on  the  surface 
pressures  is  shown.  It  is  evident  that  there  is  a  pronounced  rearward 
motion  of  the  terminating  shock  on  the  upper  surface  associated  with 
increased  blowing.  This  is  consistent  with  the  idea  that  the  latter 
cambers  the  upper  surface.  The  associated  increase  in  lift,  wave  drag, 
and  skin  friction  drag  will  be  discussed  shortly. 

In  connection  with  the  integral  thicknesses,  an  evaluation  of 
the  correctness  of  the  parabolic  solver  was  made  by  checking  their  com¬ 
patibility  with  the  integral  form  of  the  momentum  equation,  i.e.. 


(23a) 
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Fig.  5  Streamlse  development  of  boundary  layer  and  wall  jet  velocity 
profiles  over  airfoil.  Profiles  for  upper  surface  shown  In 
upper  part  of  figure;  lower  surface  values  In  lower  part  of 
figure. 
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Fig.  6  StreanMise  development  of  velocity  profiles  In  wake  for 
Moo-  0.75,  a-  1°.  Cj  »  0.055. 
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H  =  6/0 


(23b) 


T 


3U| 

%ly=0 


-  (23c) 


where  the  subscript  e  refers  to  the  edge  of  the  boundary  layer,  H  is  the 

shape  factor,  and  is  the  Mach  number  at  the  edge  of  the  boundary  layer. 

The  numerical  values  obtained  from  the  computational  solutions  satisfied 

(23)  to  within  roughly  the  truncation  error  of  the  method.  This  served 

as  a  validation  of  our  application  of  the  box  scheme  to  blown  flows  and 

also  explained  the  negative  values  of  0'(x)  on  the  blown  portions  of  the 

airfoil.  A  key  element  in  these  considerations  is  that  0  itself  is 

negative  in  these  regions,  as  has  been  previously  observed.  Thus,  if 

0  <  0  typified  by  the  upstream  portion  of  a  wall  jet,  and  an  adverse 

pressure  gradient  is  encountered,  i.e.,  f}^e  <  0,  as  in  downstream  of  the 

dx 

terminating  shock  on  the  upper  surface,  then  the  pressure  gradient  term 
P  will  be  negative.  From  (23a)  it  is  obvious  that  if 

|P|  >  Cf  . 

0'(x)  will  also  be  negative. 

Figure  4  exemplifies  a  situtation  in  which  0'(x)  as  well  as  6*'(x) 

is  predominantly  negative  over  most  of  the  blown  surface.  Since  the  skin 

friction  coefficient  can  be  influenced  by  the  blowing  coefficient  Cj, 

cases  can  be  conceived  where  the  inequality  (24)  is  reversed.  This  is 

also  possible  for  turbulent  flow.  An  an  illustration,  the  algebraic 

23  24 

edcty  viscosity  model  of  Patankar  and  Ramaprian  given  by  Eqs.  (15) 

was  used  to  compute  the  turbulent  analogue  of  the  laminar  flow  considered 

herein.  Figure  8  compares  the  chordwise  pressures  from  laminar  flow  with 

the  turbulent  case,  for  following  conditions:  M^  =  0.7,  a  »  1°,  C,  »  0.047, 

00  ^ 

Re  ■  10®  and  XsL  «  0  .19.  The  inviscid  solution  is  also  displayed.  In 
this  mildly  supercritical  situation,  it  is  evident  from  the  figure  that 
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Fig.  8  Comparison  of  laminar  and  turbulent  blown  chordwise  pressures 
with  Inviscid  theory  Moo"  0.7,  a*  1°,  Cj  •  0.047,  Re  ■  10®, 


X.,  •  0.19,  (NACA  0012  airfoil). 
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the  laminar  lifting  pressures  are  greater  than  those  predicted  by  the 
Inviscid  model.  By  contrast,  the  turbulent  values  are  less,  large 
differences  occurring  on  the  upper  surface  In  the  supersonic  zone.  Th^ 
Implication  of  this  state  of  affairs  with  respect  to  overall  forces  Is 
shown  In  Table  1. 

TABLE  1.  COMPARISON  OF  OVERALL  FORCES  ON  BLOWN  NACA  0012  AIRFOIL, 

M^  =  0.7,  a  =  1°.  Cj  =  0.047,  Re  =  10‘,  =  0.19 


Mode 

^L 

^0 

•^WAVE 

Lami nar 

0.326 

0.00716 

0.0195 

Turbulent 

0.162 

0.0129 

0.0344 

"Inviscid" 

0.287 

0.00611 

— 

To  compute  C^  »  the  nose  contribution  near  x  =  0  Is  extremely  Important. 
•^WAVE 

Although  the  singular  character  discussed  In  connection  with  Eq.  (20)  could 
be  used  for  this  purpose,  the  numerical  procedure  given  In  Ref.  29  was 
found  to  be  more  convenient.  The  trends  shown  In  Figure  8  and  Table  1  can 
be  more  easily  appreciated  by  reference  to  Figure  9  which  displays  the 
6*(x)  distributions  for  the  laminar  and  turbulent  cases.  On  comparison 
of  these  variations.  It  Is  Immediately  obvious  that  the  turbulent  flow 
through  the  attendant  larger  values  of  C^  will  at  moderate  distances  down¬ 
stream  of  the  slot  Xj|_  be  governed  by  the  Inequality 

IPj  <  .  (25) 

the  reverse  of  (24).  According  to  (23a),  e'(x)  becomes  quickly  positive, 
and  since  H  «0.5  or  less,  5*'  exhibits  a  similar  behavior*.  By  contrast 


*Downstream,  both  terms  of  the  right  hand  side  of  (23a)  become  positive. 
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Comparison  of  6*(x)'s  for  turbulent  and  laminar  flows.  Moo®  0. 
1°,  Cl  *  0.047,  Re  ®  10°,  x.,  =  0.19,  (NACA  0012  airfoil). 
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the  laminar  case  as  well  as  the  one  associated  with  Figure  4  satisfies 
(24),  and  accordingly,  6*  as  well  6*'  are  negative.  Thus,  the  turbulent 
flow  starts  with  an  excess  layer,  but  quickly  transitions  to  a  defect 
layer  with  its  larger  separation  resisting  shear  stress  at  the  wall,  in 
contradistinction  to  the  laminar  case,  which  maintains  the  excess  wall 
jet  profile  to  the  trailing  edge  with  smaller  surface  shear  stress. 

This  increased  leads  to  the  higher  Cq  shown  in  Table  1.  Because 
£*'(x)  becomes  positive  quickly  downstrelm  of  the  slot, the  more  separation 
resistant  turbulent  wall  jet  gives  rise  to  the  decambering  effect  and 
attendant  lift  loss  described  earlier.  Furthermore,  the  profile  thicken¬ 
ing  associated  6*'(x)  >  0  results  in  greater  wave  drag  than  the  laminar 
case  as  also  shown. in  Table  1. 

Returning  to  the  laminar  flow  situation,  and  the  effect  of  the 
viscous  interaction  processes  on  the  overall  force  characteristics  of 
blown  airfoils.  Figure  10  indicates  the  variations  of  wave  drag  coeffi¬ 
cient  Cn  and  friction  drag  Cr.  as  a  function  of  the  blowing  coeffi- 
‘^WAVE  “^f 

dent  Cj.  Various  flight  conditions  are  indicated,  involving  Mach 

numbers  of  0.7  and  0.75  and  incidences  a  of  1°  and  2°,  and  a  Reynolds 

number  based  on  the  airfoil  chord  of  10^  as  well  as  a  slot  location  in 

units  of  the  airfoil  chord  x^j^  =  0.19.  A  striking  and  surprising 

linearity  of  Cp.  with  C.  is  apparent  from  these  results.  This  behavior 
Uf  J 

warrants  further  stu(iy.  Comparable  variations  of  C.  appear  to  provide 

J 

considerably  lowe'"  increases  in  wave  drag  C,.  at  a  fixed  M  ,a  combina- 

“WAVE 

tion.  This  trend  which  is  believed  to  be  typical  for  other  airfoils  is 
more  clearly  demonstrated  in  Figure  11  where  the  arrows  there  and  subse¬ 
quent  plots  indicate  increasing  Cj.  From  a  practical  engineering  viewpoint, 
these  numerical  studies  corroborate  experimental  evidence  of  a  friction  or 
"scrubbing"  drag  penalty  associated  with  the  separation  delay  or  removal 
benefit  of  blowing.  The  computational  model  is  general  enough  to  provide 
cheap,  systematic  tradeoff  studies  to  optimize  a  given  airfoil.  In  this 
connection,  the  results  of  Figures  10  and  11  exemplify  another  important 
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f  F1g.  10  Effect  of  blowing  on  drag  coefficients. 
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trend,  namely,  the  substantial  increase  In  wave  drag  associated  with 
Increase  In  Mach  number  and  Incidence  at  a  given  blowing  coefficient. 

To  obtain  these  results,  a  sequence  of  blowing  coefficients  was 
Investigated  which  provided  Increasing  lengths  of  separation-free  flow 
over  the  airfoil  until  the  trailing  edge  was  reached.  Defining  the 
separation  point  x  =  at  which  the  shear  stress  vanishes  or  reverse 
flow  Is  first  encountered,  the  blowing  coefficient  was  increased  until 
the  separation  point  reached  the  trailing  edge,  x  =  1.  For  x.^p  <  1,  the 
weak  viscous  interaction  model  described  In  this  report  would  not  be 
applicable  to  describe  the  flow.  In  Ref.  6,  the  behavior  of  x^^p  with  Cj 
was  Indicated.  Further  studies  of  this  trend  were  performed  in  this 
effort  and  some  typical  results  are  shown  in  Table  2  for  a  Reynolds  number 
based  on  the  airfoil  chord  of  10®  for  various  slot  locations. 

TABLE  2.  EFFECT  OF  BLOWING  ON  SEPARATION  POINT,  M„  =  0.8,  a  =  1° 


6.2 

0.047 

0.14 

0.859 

0.796 

0.793 

8.3 

0.03 

0.14 

0.793 

0.766 

0.732 

8.1 

0.0355 

0.19 

0.84264 

0.766 

0.782 

9.1 

0.04 

0.19 

0.84264 

0.797 

0.782 

6.1 

0.0398 

0.10 

0.81890 

0.797 

0.753 

8.2 

0.0257 

0.10 

0.81890 

0.766 

0.753 

Indicated  in  the  table  for  x^^p  are  locations  of  the  x  grid  point 
closest  to  the  separation  location.  This  is  evident  in  the  computations 
when  a  divergence  of  the  Newton  iterations  employed  in  the  box  scheme 
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occurs.  This  divergence  is  associated  with  incipient  upstream  flow  or 
vanishing  shear  stress.  Also  tabulated  are  the  locations  of  the 
upstream  (low  pressure  side)  of  the  shock  layer  (which  was  spread  ih  ^ 
the  inviscid  zero^^  iterate  over  two  to  three  points  in  the  SOR  grid),  and 
Xj^,  the  X  location  of  the  last  box  scheme  point  upstream  of  x^^p.  From 
the  results,  it  is  evident  that  separation  occurred  close  to  the  foot  of 
the  shock  in  these  cases.  The  superficially  puzzling  invariance  of  x^^p 
to  moderate  changes  in  C.  for  Cases  8.1,  9.1,  6.1  and  8.2  can  be  explained 
on  the  basis  of  the  relative  coarseness  of  the  box  scheme  grid  near  the 
shock  foot.  It  is  evident  from  Table  2  that  with  the  exception  of 
Case  8.1  where  the  shock  was  upstream  of  the  last  box  scheme  grid  point 
prior  to  separation,  the  shock  was  between  this  point  and  x^^p.  With  the 
relative  coarseness  of  the  box  scheme  grid  in  this  location,  it  is  con¬ 
ceivable  that  identical  x„„„  values  could  be  obtained  in  the  shock-induced 

sep 

adverse  pressure  gradient  in  this  location. 

Defining  C’J  as  the  blowing  coefficient  at  which  x  becomes 
J  sep 

unity  as  C.  from  below,  values  of  C.  »  C^  ("overblowing")  may  be 

«J  J  J  V 

inefficient  from  an  L/D,  Cr,  and  frictional  or  scrubbing  drag  viewpoint, 

‘^WAVE 

despite  supercirculation  gains.  In  addition  to  Figure  10  and  11,  the 
importance  of  frictional  drag  in  relation  to  wave  drag  is  shown  in  Figure  12. 
Although  large  increases  in  the  friction  drag  were  encountered  with  increases 
in  C.,  the  percentage  of  this  mode  to  the  total  is  not  strongly  influenced 

J 

by  the  magnitude  of  the  blowing  coefficient.  It  is  evident  that  the  lower 
Mach  number- incidence  combinations  lead  to  frictional  drags  of  the  order 
of  70  to  75%  of  the  total,  whereas  for  a  modest  increase  in  Mach  number, 
it  reduces  to  about  half  of  the  total.  To  our  knowledge,  these  are  the 
first  quantifications  of  such  trends  from  a  computational  model. 

The  variations  of  lift  coefficient  associated  with  overblowing  are 
shown  in  Figure  13.  Relatively  small  increases  in  lift  are  exhibited  for 
appreciable  increases  in  blowing.  Such  increases  can  lead  to  degradations 
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Cj  -  BLOWING  COEFFICIENT 


Fig.  12  Fraction  of  total  drag  associated  with  various  blowing  rates 
and  flight  conditions. 
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0.054  0.058 


Fig.  13  Effect  of  blowing  on  lift 
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in  aerodynamic  efficiency  (L/D)  as  shown  in  Figures  14  and  15  for  these 
blowing  enhancements.  Further  elucidation  of  these  trends  is  shown  in 
Figure  16,  where  the  L/D  changes  associated  with  lift  gains  are  shown 
(setting  aside  the  issue  of  trim  drag).  For  the  lower  Mach  number  cases, 
within  the  Cj  range  assumed  (which  is  somewhat  greater  than  C^),  increases 
in  lift  give  L/D  loss.  At  the  higher  Mach  number,  an  optimum  is  shown 
for  maximum  L/D.  It  is  interesting  to  note  that  this  is  obtained  at  the 
smallest  C-  studied,  i.e.,  0.0355.  Although  the  achievement  of  high  lift 

J 

is  necessary  for  good  turn  performance,  the  significance  of  the  L/D  loss 
is  important  from  the  viewpoint  of  maintenance  of  constant  altitude,  speed, 
or  specific  power  factor  during  maneuver  scenarios. 

In  addition  to  the  foregoing  aspects,  the  viscous  interaction 
model  was  used  to  study  the  effect  of  Reynolds  number  and  slot  location  on 
aerodynamic  characteristics  as  shown  in  Table  3*. 


TABLE  3.  EFFECT  OF  REYNOLDS  NUMBER  Re  AND  SLOT  LOCATION  ON 
AERODYNAMIC  CHARACTERISTICS,  M  =  0.75,  a  =  1°, 

C.  =  0.0371 

J 


Re 

^SL 

\ 

Cl 

10® 

0.10 

0.0142 

0.0177 

0.531 

10^ 

0.10 

0.0112 

0.0144 

0.360 

10’ 

0.25 

0.0153 

0.0152 

0.383 

With  the  increased  Reynolds  number,  the  friction  drag  and  lift 
coefficient  are  reduced.  This  is  in  accord  with  the  transpiration 
boundary  condition  modification  effects  previously  described.  The  reduced 
lift  is  associated  with  a  reduction  of  viscous  camber  and  supercriticality 
at  the  higher  Reynolds  number.  Associated  diminished  profile  thinning 

♦For  X5j_*0.10,  the  "inviscid"  results  gave  Cq^»  0.00863  and  C|_  »  0.322, 
consistent  with  the  Re  trend  shown  in  Table  3. 
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reduces  the  wave  drag.  Table  3  also  shows  reduced  friction  drag  associated 
with  downstream  movement  of  the  slot  for  the  same  blowing  coefficient. 

This  Is  related  to  a  smaller  coverage  of  the  enhanced  frictional  dlsslga- 
tlon  of  the  wall  jet  from  the  slot  blowing.  Since  the  net  sink  strength 
Is  less  In  proportion  to  the  region  covered,  the  effective  thickness  is 
greater,  leading  to  increased  wave  drag  and  reduced  lift  with  rearward 
slot  movement. 

From  the  model  and  the  studies  previously  described,  certain 
considerations  emerge  relevant  to  the  design  of  optimum  blown  transonic 
airfoils.  Important  ones,  we  believe,  are  the  following: 

1.  In  addition  to  the  possibility  of  high  friction  drag 

occurring  with  excessive  overblowing,  additional  wave  drag 
arises  due  to  a  cambering  effect  from  the  downward  momentum 
imparted  to  the  flow  from  the  wall  jets.  To  deal  with  this 
problem.  Inverse  methods  such  as  those  developed  In  Ref.  19 
generalized  to  the  viscous  Interaction  formulation  given 
herein  could  lead  to  Improved  aerodynamic  efficiencies 
associated  with  lift  augmentation  with  blowing.  The  Inverse 
procedure  would  be  utilized  to  eliminate  or  suppress  shock 
waves  In  the  flow.  In  addition  to  the  attendant  wave  drag 
reduction,  there  would  also  be  an  accompanying  alleviation 
of  the  adverse  pressure  gradient  creating  separated  flow. 

This  modification  would  reduce  the  C.  required  to  achieve 
attached  flow.  In  turn  easing  engine,  structural,  and  heat 
transfer  penalties  as  well  as  offsetting  internal  flow 
losses.  Conversely,  although  Increased  blowing  delays 
separation,  there  is  a  secondary  counteractive  effect 
Involving  the  intensification  of  the  shock  (due  to  the 
induced  camber),  augmenting  the  adverse  pressure  gradient. 
This  will  degrade  the  effectiveness  of  the  blowing  to  some 
extent. 
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2.  The  possibility  of  modifying  transition  and  achieving 

laminar  flow  through  tangential  blowing  Is  of  some  Interest. 
One  conceivable  mechanism  Is  through  alteration  of  the  ^ 
velocity  profiles.  For  unblown  cases,  Rayleigh's  Inviscid 
stability  theory  (see  Ref.  30  for  an  account)  Indicates  that 
a  necessary  condition  for  amplification  of  disturbances  Is 
the  occurrence  of  a  point  of  Inflection  In  the  velocity  pro¬ 
file.  In  the  case  of  monotonic  velocity  profiles,  Rayleigh's 
Inflection  point  theorem  can  be  derived  directly  from  the 
Sturm-Llouville  problem  related  to  a  plane  wave  stability 
analysis.  Based  on  this  thinking,  pressure  gradient  param¬ 
eters  and/or  the  shape  factor  H  defined  by  the  ratio  of  the 
displacement  to  momentum  thickness  can  be  correlated  against 
the  critical  Reynolds  number  R^  for  the  Inception  of  Insta¬ 
bility.  In  this  context.  It  has  been  demonstrated  theoret¬ 
ically  and  experimentally  that  favorable  pressure  gradients 
greatly  enhance  stability  and  unfavorable  ones  detract  from 
It.  Moreover,  decreasing  H  toward  2  from  about  3.6  has  been 
shown  to  substantially  Increase  R^. 


Application  of  Item  2  to  tangential  blowing  and  wall  jets  has 
not  been  made  to  our  knowledge,  and  leads  to  some  difficulties.  Super¬ 
ficially,  It  might  be  anticipated  that  blowing  could  be  used  to  control 
the  shape  factor  and  through  generalizations  of  the  considerations  asso¬ 
ciated  with  boundary  layers,  augment  stability,  with  the  ultimate  appli¬ 
cation  of  lamlnarizatlon  of  an  airfoil.  Indications  of  the  alterations 
of  the  curvature  of  the  velocity  profiles  and  changes  In  the  shape  factor 
associated  with  the  blown  airfoil  considered  In  this  Investigation  are 
given  In  Figures  17  and  18.  It  Is  Interesting  to  note  that  reflecting 
the  excess  velocity  character  of  the  wall  jets  on  these  airfoils,  the 
shape  factors  shown  In  Figure  18  are  nowhere  In  the  range  considered  In 
the  standard  stability  literature.  In  violation  of  the  conditions  of 
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Fig.  17  Curvature  of  velocity  profile  at  wall  for  varying  blowing 
Mod-  0.75, «»  1  ,  (upper  surface). 
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Rayleigh's  theorem,  the  velocity  profiles  of  wall  jets  are  furthermore 
non-monotone  in  character.  An  additional  set  of  practical  complications 
deals  with  density  variations  in  the  hot  engine  gas  used  for  the  tangen¬ 
tial  blowing  as  well  as  upstream  disturbances  that  could  act  as  trips 
upon  injection  into  the  boundary  layers.  In  spite  of  these  difficulties, 
it  is  felt  that  this  area  is  a  fruitful  one  for  future  theoretical  and 
experimental  research.  Regarding  the  former,  a  useful  point  of  departure 
would  be  investigation  of  the  stability  properties  of  simple  wall  jets 
according  to  a  generalization  of  the  thinking  leading  to  Rayleigh's 
theorem. 
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5.0  CONCLUSIONS 

A  viscous  interaction  model  has  been  developed  for  the  treatment 
of  tangentially  blown  transonic  airfoils.  The  coupling  procedure  employs 
a  small  disturbance  approximation  for  the  nearly  inviscid  outer  flow  and 
interacts  this  region  with  the  boundary  layers,  wall  jets  and  propulsive 
wake.  The  box  scheme  has  been  generalized  to  treat  the  latter,  giving 
rapid  convergence  of  the  viscid-inviscid  iterative  sequence.  Observations 
based  on  the  parametric  studies  conducted  indicate  that  for  laminar  flow: 

1.  The  skin  friction  drag  appears  to  increase  nearly  linearly  with 
blowing  coefficient  Cj  for  a  fixed  Mach  number- incidence 
combination  giving  a  substantial  "scrubbing"  drag  penalty 

for  overblowing  beyond  the  Cj  =  C^J  necessary  to  achieve 
attached  flow. 

2.  Wave  drag  increases  are  much  smaller  than  those  for  skin 
friction  over  the  same  range  of  increase  of  Cj. 

3.  The  viscous  interactions  move  the  terminating  shock  of  the 
supersonic  region  rearward  with  respect  to  its  inviscid 
position.  There  is  an  attendant  increase  in  lift  and 
supercriticality  with  the  viscous  interactions. 

4.  The  computational  model  can  quantify  the  rearward  movement 
of  the  separation  point  with  increases  in  Cj  and  thereby 
be  used  to  size  the  thrust  bleed  requirements  for  mainte¬ 
nance  of  attached  flow.  Other  trade-offs  are  also  easily 
accessible  using  the  simulation. 

5.  Although  large  increases  in  friction  drag  were  encountered 
with  moderate  increases  in  Cj  at  a  given  Mach  number  and 
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incidence,  the  percentage  of  the  friction  drag  to  the 
total  drag  is  not  strongly  influenced  by  Cj  over  the 
range  considered.  . 

6.  For  a  fixed  Cj,  the  percentage  of  friction  drag  to  the  total 
drag  tends  to  decrease  at  higher  Mach  number-incidence 
combinations,  presumably  reflecting  the  attendant  dispropor¬ 
tionate  increases  in  wave  drag. 

7.  The  supercirculation  gains  associated  with  overblowing 

(Cj  >  CJ)  do  not  provide  large  increases  in  the  aerodynamic 
efficiency  (L/D)  due  to  the  scrubbing  drag  referred  to  in 
Item  1. 

8.  Increased  Reynolds  number  leads  to  decreases  in  the  lift 
augmentation  associated  with  viscous  interaction  effects  as 
well  as  reduced  friction  and  wave  drag. 

9.  Downstream  slot  movement  increases  the  wave  drag,  as  well 
as  reducing  the  skin  friction  drag  and  lift.  This  is 
particularly  important  from  a  design  viewpoint  with  respect 
to  minimization  of  the  scrubbing  drag. 

10.  Studies  based  on  an  algebraic  turbulence  model  employed  in 
the  viscous  interaction  algorithm  suggest  that  the  lift 
decreases  from  values  predicted  by  the  inviscid  model  in 
contrast  to  the  increases  in  the  laminar  case.  There  is 
also  an  attendant  increase  in  frictional  and  wave  drags. 

This  behavior  can  be  explained  from  the  presence  of 
greater  friction,  which  reverses  the  laminar  negative 
streamwise  gradient  of  displacement  thickness  downstream 
of  the  injection  slot,  decambering  the  airfoil  as  in  the 
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case  of  a  conventional  boundary  layer  and  in  contrast  to 
the  cambering  effect  of  laminar  blowing.  The  reversal, 
which  is  observed  from  the  computational  solutions  can  be. 
explained  in  terms  of  the  Karman  integral  momentum 
equation.  Thus,  the  enhanced  resistance  to  separation  within 
the  turbulent  framework  is  accompanied  by  a  degradation  of 
aerodynamic  characteristics  as  compared  to  laminar  flow. 
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6.0  RECOMMENDATIONS 

On  the  basis  of  the  foregoing  discussion,  it  is  evident  that  4. 
number  of  important  areas  should  be  addressed  to  understand  and  exploit 
tangential  blowing  procedures  for  the  improvement  of  the  aerodynamic 
performance  of  wings.  From  a  modeling  viewpoint,  the  scheme  discussed 
applies  certain  approximations  at  the  slot,  trailing  edge,  and  far  field. 

In  the  vicinity  of  the  slot,  the  boundary  layer  assumptions 
break  down  and  wave  interactions  can  play  a  key  role  in  providing  the 
appropriate  initial  conditions  for  the  ensuing  downstream  boundary  layer 
development.  In  this  vicinity,  waves  may  reflect  or  transmit  through 
the  mixing  layer  modifying  the  external  flow  and  the  wall  jet.  The 
asymptotic  procedures  formulated  in  Ref.  6  should  be  further  generalized 
and  computationally  applied  to  treat  these  interactions.  This  analysis 
could  be  useful  in  determining  the  slot  location  optimizing  separation 
prevention  and  aerodynamic  characteristics. 

Near  the  trailing  edge,  local  solutions  which  are  essentially 
a  generalization  of  existing  multiple-deck  methods  for  blowing  should  be 
developed.  The  implications  of  this  refinement  with  respect  to  the 
Kutta  condition  and  lift  augmentation  should  be  assessed.  In  spite  of 
recent  work  on  the  near  wake,  there  appears  to  have  been  no  treatments 
of  cases  in  which  the  shock  is  near  the  trailing  edge  or  the  latter  is 
close  to  critical  conditions  even  for  the  unblown  case.  This  situation 
occurs  frequently  in  practice,  and  a  detailed  analysis  of  it  generalized 
to  blowing  would  have  great  impact. 

The  analysis  in  this  report  utilizes  a  far  field  based  on 
circulation  theorems  for  viscous  two-dimensional  flow  and  heuristic 
considerations.  A  systematic  asymptotic  treatment  involving  matching, 
highlighting  the  role  of  viscous  effects  on  the  vorticity  in  the  wake 
as  well  as  the  development  of  the  supercirculation  field  with  blowing 
would  be  an  extremely  useful  analysis.  It  would  provide  further  insight 
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into  the  mechanisms  associated  with  thrust  recovery  as  well  as  Illuminate 
the  structure  of  blown  flow  fields. 

From  a  practical  viewpoint,  the  simultaneous  use  of  an  Inverse* 
method  with  the  viscous-interaction  model  could  substantially  offset 
the  wave  drag  penalty  associated  with  blowing  strong  enough  to  prevent 
separation.  Further  effort  along  these  lines  Is  therefore  recommended. 

In  addition  to  the  foregoing,  application  of  blowing  an  three- 
dimensional  configurations  leads  to  other  Important  Issues.  Some  of 
these  are  related  to  further  compromises  required  to  deal  with  vortex 
drag  In  addition  to  wave  and  friction  forms,  as  well  as  the  possibilities 
of  exploitation  of  favorable  Interference  effects  interacting  with  those 
from  blowing  to  optimize  aerodynamic  performance.  Examples  of  the  latter 
are  blown  canard-wing  arrangements.  In  whicf^  :he  blowing  can  offset 
adverse  pressure  gradients  Induced  by  the  canard  Induction  on  the  wing 
transonically  and  yet  preserve  low  speed  stall  alleviation  through  vortex 
Interactions. 
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